Numerical investigation of a typhoid disease model in fuzzy environment

Salmonella Typhi, a bacteria, is responsible for typhoid fever, a potentially dangerous infection. Typhoid fever affects a large number of people each year, estimated to be between 11 and 20 million, resulting in a high mortality rate of 128,000 to 161,000 deaths. This research investigates a robust numerical analytic strategy for typhoid fever that takes infection protection into consideration and incorporates fuzzy parameters. The use of fuzzy parameters acknowledges the variation in parameter values among individuals in the population, which leads to uncertainties. Because of their diverse histories, different age groups within this community may exhibit distinct customs, habits, and levels of resistance. Fuzzy theory appears as the most appropriate instrument for dealing with these uncertainty. With this in mind, a model of typhoid fever featuring fuzzy parameters is thoroughly examined. Two numerical techniques are developed within a fuzzy framework to address this model. We employ the non-standard finite difference (NSFD) scheme, which ensures the preservation of essential properties like dynamic consistency and positivity. Additionally, we conduct numerical simulations to illustrate the practical applicability of the developed technique. In contrast to many classical methods commonly found in the literature, the proposed approach exhibits unconditional convergence, solidifying its status as a dependable tool for investigating the dynamics of typhoid disease.

including infection resistance, Nthiiri et al. developed a mathematical model.Using the next-generation matrix technique, the model's stable states are identified and the reproduction number is calculated.The model's stability is analyzed to identify the factors that contribute to the disease's spread within a specific community 5 .The dynamics of the typhoid fever model were investigated, as well as the existence and uniqueness of the solution, by Peter et al.For the model, stability analysis is also carried out 6 .Bakach et al. reviewed some mathematical models of typhoid 7 .Karunditu et al. formulated a mathematical model of typhoid fever incorporating unprotected humans.The local and global stability of equilibrium points is also studied 8 .A mathematical model for the transmission of typhoid was developed by Nyaberi and Musaili, and it examines the effects of treatment on the dynamics of the illness 9 .Birger et al. studied mathematical models of typhoid transmission by considering FQNS and multidrug resistance separately.The effect of vaccination was predicted on the basis of forecasts of vaccine coverage 10 .A mathematical model for typhoid fever spread in a population is formulated.The equilibrium points of the model and their stabilities are investigated 11 .By utilizing several optimal control strategies, Wameko et al. established a mathematical model to look into the dynamics of typhoid disease.Typhoid disinformation among the population is reduced when the three control techniques are quickly implemented, as demonstrated 12 .
Due to population-wide variations in susceptibility, exposure, infectivity, and recovery, the parameters employed in epidemic models are imprecise.If different age groups, population groups, and resistance patterns are taken into account, differences may result.To take into consideration these varying degrees of persons, more realistic models are required.Mishra et al. claim that due to the high level of uncertainty, epidemic systems, particularly those involving infectious diseases, require a new approach 13 .Fuzzy sets and fuzzy logic have been extensively utilized to tackle real-world problems across diverse domains, encompassing medicine, engineering, economics, and numerous other fields where human decision-making plays a pivotal role in assessment and logical reasoning [14][15][16][17][18] , just to mention a few.Moreover, scholars have harnessed this theoretical framework in epidemiology as well.Incorporating fuzzy theory and treating the transmission coefficient as a fuzzy set, Barros et al. suggested a SI model 19 .Fuzzy logic was used by Ortega et al. to predict issues with infectious disease epidemiology.A rabies model in dogs with incomplete vaccinations was discussed 20 .
Mondal et al. developed an SIS model for investigating the plague using the fuzzy set theory 21 .Das and Pal developed a SIR model and studied it mathematically and numerically 22 .Sadhukhan et al. conducted research on harvesting optimization in a food chain model in a fuzzy environment 23 .To capture the dynamics of coronavirus illness, Li et al. developed a fuzzy SEIR model supplemented by confidence index theory 24 .Abdy et al. presented an SIR model that incorporated fuzzy parameters to depict the dynamics of COVID-19 25 .Furthermore, Allehiany et al. explored a fuzzy SIR model employing Euler, RK-4, and NSFD methods 26 .The NSFD approach, first described by Micken 27 , has been used by a number of researchers for solving systems of differential equations [28][29][30] , to name a few.Adak and Jana investigated an SIS model involving treatment control with the utilization of fuzzy numbers 31 .
The existing mathematical models of typhoid are insufficient for the advancement of fuzzy numerical and mathematical procedures.We investigated a typhoid model with fuzzy parameters with this in mind.We can cope with the difficulties of uncertainty quantification in mathematical disease modeling by using fuzzy theory.As a result, the introduction of fuzzy parameters aids in our ability to comprehend the dynamics of typhoid transmission.Even the biological factors employed in mathematical models are not always constant because each community changes as the environment change.The majority of the issues linked to the rise in the earth's average temperature are caused by global warming.The rate at which the virus spreads throughout the population is also impacted by temperature changes.Fuzzy mathematical models are more insightful than crisp models in this regard.The creation, use, and analysis of first order explicit numerical techniques in fuzzy non-standard finite difference situations are novel aspects of the current work.
The rest of this study is organized as follows: The formulation of the fuzzy model is discussed in the following section.Within the same section, we discuss equilibrium analysis, stability analysis, and the fuzzy basic reproduction number.Following that, in the numerical modeling section, we elaborate on the creation of the forward Euler scheme and NSFD schemes for the examined model.This section also includes an assessment of the NSFD scheme's stability and consistency.The next part displays numerical simulations involving the developed techniques.Finally, in the final section, we summarize concluding observations and outline future research directions.

Typhoid fever model with fuzzy parameters
We examined the mathematical model previously discussed by Arif et al. 2 Here, T 1 represents the proportion of humans who are protected, T 2 represents the proportion of humans who are susceptible, T 3 signifies the fraction of humans who are currently infected, and T 4 indicates the fraction of (1) humans receiving treatment.The rate at which humans receiving treatment transitions from the fraction of infected humans is denoted as β .The variable α represents the rate at which individuals enroll in the protected human compartment against typhoid, while (1 − α) denotes the rate of individuals who remain susceptible to the virus, δ signifies the rate at which individuals experience a transient phase due to typhoid fever, and θ reflects the per capita rate at which susceptible individuals contract typhoid fever infection.Lastly, µ represents the natural rate of human death and birth.The depiction of the communication dynamics in the typhoid fever model is presented in Fig. 1.
The fuzzy representation of the model mentioned above can be expressed as follows: We assume that susceptible humans contract typhoid fever infection at a per capita rate denoted by θ(�) and the fraction of treated humans stemming from infected individuals, β(�) , are represented as fuzzy numbers, contingent upon the individual's virus load.The variable θ(�) can be defined as 19 The θ(�) reaches its peak when is at its maximum value, and it becomes insignificant when is at its minimum.min represents the minimum virus load required for disease transmission, and disease transmission is at its highest when equals M , reaching a value of 1.Similarly, β(�) can be defined as 32 where β 0 > 0 is the minimum treatment rate.

Equilibrium analysis
The analyzed model possesses a disease-free equilibrium point (DFE) and two endemic equilibrium points (EE).
Case 1.If � < � min and β(�) > 0 then θ(�) = 0 and we obtain, In this scenario, the human population remains free from typhoid, and this state is referred to as the DFE.From a biological perspective, typhoid disease is considered eradicated when the disease concentration within the population falls below the minimum threshold required for its sustained existence. (5)

Stability analysis
To check the stability, let us assume the following system: and The Jacobian of the system ( 11)-( 14) can be represented as: The Jacobian at the DFE is Vol.:(0123456789)The local asymptotic stability of the steady-state is confirmed when the absolute eigenvalues of the Jacobian matrix mentioned earlier having negative real parts.Analyzing the Jacobian matrix, we find that the eigenvalues are as follows: 1 = −(γ + µ) , 2 = 4 = −µ , and 3 = −(δ + β(�) + µ) .The fact that all of these eigenvalues having negative real parts validates the desired result.

Fuzzy basic reproductive number R f t
Using the next-generation matrix technique, we calculated the reproductive number, denoted as R t .
Verma et al. 33 examined different scenarios by manipulating the parameters and, in each case, determined the reproduction number to assess whether the virus spread was effectively managed.Now, R t , being a function of the typhoid disease, can be analyzed as follows 34 : . The disease-dependent function R t (�) correlates positively with the disease parameter , and its definition includes a fuzzy vari- able.As a result, the expected value of R t (�) is well defined, and its representation can be written as a triangular fuzzy number, as follows: Now, R f t can be found as and therefore,

Forward Euler scheme
The Forward Euler scheme is a well-known explicit first-order numerical approach for solving ordinary differential equations.It is computationally efficient and provides a rapid estimation of the behavior of solutions over time.

Non-standard finite difference (NSFD) scheme
The NSFD scheme is a class of numerical methods for approximating solutions to differential equations.These approaches differ from traditional finite difference methods in their approach to discretizing the domain and approximating derivatives.It has the potential to improve accuracy in the solution of differential equations.The NSFD numerical model is formulated based on the NSFD theory introduced by Mickens 35 .

Stability of the NSFD scheme
To study the stability analysis, let us assume the following system: The Jacobian matrix corresponding to the system (67)-( 70) is Jacobian at the DFE is Eigenvalues of the above Jacobian matrix are 1 = Because all eigenvalues are less than one, this validates the desired outcome that the NSFD scheme is stable at the DFE 36 .NSFD schemes do not constitute a singular category of numerical methods.Numerous researchers have expanded upon the Mickens theory.For instance, Gurski 37 introduced a straightforward mathematical framework for NSFD schemes, specifically tailored for small systems of nonlinear differential equations.This approach leverages conventional techniques used in approximating differential equations, including the incorporation of artificial viscosity and the implementation of a predictor-corrector scheme.In their work, they examined both the NSFD scheme proposed by Mickens and the one developed by Erdogan and Ozis for first-order equations. (61) (62) (64) (66)

Numerical simulations
In this section, we present the simulation results of the performance of Euler method illustrated in Fig. 2a-f, and the simulation results obtained using the NSFD method illustrated in Fig. 3a-f.The graphical representations of the performance of Euler method are illustrated at various step sizes in Fig. 2a-f.Initially, at a small step size h = 1 , the method exhibits stability, positivity, and convergence.However, as the step size is slightly increased to h = 5 , the method begins to oscillate and generates negative values in all three cases.In models like this, negative values are not meaningful since all compartments represent populations, Figure 2. The portions of sub populations using Euler method at different step sizes.and negativity is not feasible.Consequently, it can be deduced that this method is not a dependable tool for describing such models.In Fig. 3a-f, we present the simulation results obtained using the NSFD method.This time, the method demonstrates stability, both at smaller and larger step sizes, and an increase in the step size does not adversely affect its convergence and positivity, which are crucial characteristics for modeling disease dynamics.The method consistently displays stable and convergent behavior across all three cases, underscoring its reliability for studying disease dynamics in such conditions.It is worth noting that many classical numerical schemes and their fuzzy counterparts tend to lose their convergence and positivity and struggle to maintain stability as the step size increases.In contrast, the proposed method proves to be more efficient and capable of addressing such issues.We provide the convergence analysis for proposed NSFD scheme analytically and then verify it numerically in numerical simulations section.It can be noted here that the proposed NSFD scheme gives unconditional convergence and remains consistent with the continuous dynamical system.Figure 3a-f show the convergence of NSFD scheme to true equilibria, retaining all the essential features of continuous model unlike the Euler scheme which fails to do so (see Fig. 2a-f).

Conclusions
In this study, we have considered a typhoid model with fuzzy parameters.We assumed that the infection does not transmit equally among the individual of the populations.Similarly, the treatment rate is also not the same for each individual.As a function of the virus concentrations, we treated the typhoid transmission rate θ(�) and the treatment rate β(�) as fuzzy variables.In deterministic models, these parameters are fixed and independent of the viral load.As a result, it may be said that the fuzzy model is more adaptable and balanced than the crisp system.Fuzzy theory is used to address uncertainty quantification difficulties in mathematical disease modeling.We examined it for various virus loads because fuzzy variables are functions of virus load that depend on virus levels.In light of this, we addressed the studied model's fuzzy equilibrium points while taking the population's virus levels into account.We established that the disease-free equilibrium point is reached if the virus concentration is lower than the minimal concentration necessary for disease transmission in the population.When the population's viral levels exceeded the bare minimum needed for disease transmission, we reached the endemic equilibrium points.For various viral concentrations, the basic reproduction rate is examined.We employed the expected value of a fuzzy number to ascertain the fuzzy basic reproduction number.Two numerical schemes are developed for the approximate solution of the studied model.The developed schemes are analyzed for different amounts of virus.The suggested numerical algorithms must maintain the positive nature of the solutions of such dynamic population models.The Euler method preserved this for only small values of the step sizes and generated negative values by increasing the step size.On the other hand, the NSFD preserved this for all large values of the step size too for different amounts of virus.Additionally, the convergence and consistency of the NSFD scheme are analyzed, demonstrating that the suggested approach is unconditionally convergent and consistent of order 1.The creation, application, and analysis of a non-standard finite difference technique for the numerical analysis of a typhoid illness model with fuzzy parameters are the main foci of the current work.Future developments may include fuzzy stochastic, fuzzy delayed, or fuzzy fractional dynamic senses.Age-structured fuzzy epidemic models could potentially be modeled using the NSFD modeling theory.This research focuses mostly on using triangular fuzzy numbers as membership functions.In the future, we may investigate the use of various fuzzy numbers as potential membership functions, such as trapezoidal, pentagonal, hexagonal, and so on.

Figure 3 .
Figure 3.The portions of sub populations using NSFD method at different step sizes.